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are discussed. The finite volume corrections to the probability distributions and particle density 
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I. INTRODUCTION 

Statistical models have been shown to be very successful in the description of particle production in heavy ion 
collisions [3]. Such models are usually constructed from the partition function of a quasi non-interacting gas composed 
of all known hadrons including hadronic resonances. The contribution of resonances is an effective approach to 
reproduce strong interactions in the hadronic medium (2]. The hadron resonance gas model has also been shown to 
be consistent with lattice QCD thermodynamics restricted to the confined, hadronic, phase Mi 

Considering a hadronic system using statistical models, it is essential to implement constraints related to internal 
symmetries Jij, ij. We consider here an abelian symmetry corresponding to one conserved charge. In the statistical 
system, conservation laws can be implemented in the canonical (C) or grand canonical (GC) ensembles. In the 
following we consider only ultrarelativistic systems as encountered e.g. in high energy heavy ion collisions. In such a 
case particle numbers are not conserved, thus there is no associated chemical potentials to the particle numbers. The 
only chemical potentials included here are those related to conserved charges. In the grand canonical ensemble (GC) 
the quantum numbers have fixed average values which are determined by the corresponding chemical potentials, on 
the other hand, in the canonical ensemble, quantum numbers have fixed values. This leads to an essential difference 
in the volume dependence of observables in the GC and C formulations musoa In the limit when V —> oo some 
ratios of extensive quantities converge to different values in the GC and C ensembles 0 12j. The thermodynamic 
limit (T-limit) is realized if the volume V —> oo while the charge density q and the particle densities rii remain finite 
in the C system while their thermal average values (q) and (m) in the GC system are kept fixed. It is thus clear 
that the equivalence of both descriptions in the thermodynamic limit can be strictly established only for intensive 
observables m. The equivalence of the GC and C descriptions in the thermodynamic limit has been established at 
the most basic, probability density level. It has been shown recently [^j that different particle probability densities 
calculated from the C and GC ensembles coincide in the thermodynamic limit. 

In some physical situations, however, the knowledge of T-limits of probability density functions is not sufficient. 
There is a broad class of physical quantities which are finite in the thermodynamic limit but are still different for 
different statistical ensembles. We call henceforth such variables semi-intensive quantities. These variables have 
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finite T-limits, like intensive variables, but those limits are governed by finite volume corrections to probability 
distributions for particle densities. The properties of semi-intensive quantities in the T-limit are entirely determined 
by the 0(V -1 ), subleading corrections to the probability density functions. These contributions are also important 
if one studies finite volume corrections to thermodynamic observables. Such a situation appears when comparing the 
statistical model with lattice gauge theory results obtained on a small lattice. The statistical model description of 
finite volume effects e.g. the charge susceptibilities calculated recently on the lattice HQ require knowledge of the 
finite volume corrections to the probability densities. 

The phenomenologically relevant example of a semi-intensive quantity is the scaled variance 


_ AN 2 _ (N 2 ) - ( N) 2 

UJn -Jn)- W) 


( 1 ) 


where N is, e.g., the number of charged particles. The u>n is an experimental observable in heavy ion collisions and 
measures the relative charged particles fluctuation in a system Hi- 

it was recently pointed out [llj that the scaled variance m has different thermodynamic limits in the GC and C 
ensembles. In fact, by the substitution (N k ) = V k (n k ), the ojn is related to the scaled variance for charged particle 
densities to n as 


AN 2 V 2 An 2 rr 

UJN ~ ~{Ny ~ V(n) = Vw ' 


( 2 ) 


In the T-limit the scaled variance u> n vanishes both in the GC and in the C ensembles. We show that this can be 
formally expressed as 

uj n = ^n + o(y- 2 ). (3) 

Consequently the T-limit of the scaled variance u>n is sensitive to the next to leading order (NLO) corrections 0( Id -1 ). 
The subleading term 7Z is specific to a given ensemble, that is why the T-limit of a; at is ensemble dependent. 

In this paper we calculate the NLO corrections to different probability distribution functions in the C and the GC 
ensembles. These subleading corrections will be explicitly determined in the hadron resonance gas model constrained 
by charge conservation. We will show that in the T-limit the probability density distributions do not exist as regular 
but appear as generalized functions. We also calculate the T-limit properties of the related particle density moments. 
Finally, as one of the applications, we introduce a class of semi-intensive variables in the canonical and the grand 
canonical ensemble and discuss their behavior in the vicinity to the thermodynamic limit. 


II. DISTRIBUTIONS AND MOMENTS IN FINITE SYSTEMS 


The properties of charged particle probability distributions in their approach to the thermodynamic limit will be 
discussed in the context of the statistical model of a non-interacting gas constrained by the conservation of the 
abelian charge Q. The thermodynamic system of volume V and temperature T is considered to be composed of 
charged particles and their antiparticles carrying charge ±1 respectively. The requirement of charge conservation in 
the system is imposed on the grand canonical or canonical level. 

The partition function of the above C and GC statistical system is found to be 


Z$(V,T)= Tr Q e-P* = I Q (2z) , 

Z gc (V,T ) = Tre - ^ - ** 45 = exp (2zcosh^) . 

where z is the sum over all one-particle partition functions 


z(T) 


(2tt)3 4 ^ 9i 


d 3 p e -/3 VP 2 + m2 i 


m 2 K 2 (^) = Vz 0 (T ), 

i 


(4a) 

(4b) 


( 5 ) 


and gi is the spin degeneracy factor. The sum is taken over all charged particles and resonances of mass mi carrying 
the charge ±1. The functions Iq and I \ 2 are modified Bessel functions. The chemical potential p determines the 
average charge in the GC ensemble. 

Semi-intensive variables are constructed from different particle moments ( N k ) and their volume dependence is 
obtained from the corresponding behavior of the particle number probability distribution. 
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In the C ensemble of a system of volume V and total charge Q the probability distribution P%(N, V) to have N 
negatively and N + Q positively charged particles is obtained mm from the partition function (14all as 


y 2N-\-Q 


pQ(iV,y ) N\(N + Q)\ Iq(2z) ' 


( 6 ) 


On the other hand in the GC ensemble with volume V and average charge (Q) the probability distribution 
P(q\ (IV, Q, V) to find a system with a given charge Q and a given number of negatively charged particles N is 
expressed m as the product 

VfQ) (N, Q, V) = (TV, V) V${Q, V ), (7) 


of the canonical particle number distribution Vq(N, V ) from Eq. © and the grand canonical probability distribution 


'P?q ) (Q,V) = Iq(2z) 


('Q) + y/(Q) 2 +4z 2 
2z 


e ^ v /(Q) 2 + 4 ^ 


( 8 ) 


to find the total charge Q in the system with the average charge (Q). 

With the knowledge of the probability distributions from Eqs. © and © the thermal average of particle moments 
in the GC and C ensemble are obtained from 

OO OO 

(N k )f Q C } = E E N k V? Q c } (N,Q,V), 

N=0Q=-N 
oo 

(N% = E^qW^)' 

N=0 

In the following section we will discuss the generalizations of the above results for the probability distributions and 
particle moments that are required to analyze the approach to the thermodynamic limit. 


(9a) 

(9b) 


III. DISTRIBUTIONS AND MOMENTS IN INFINITE SYSTEMS 

In the previous section we have summarized how to relate particle moments with probability distributions in a 
system that is constrained by charge conservation. These results, however, are only valid for finite systems far away 
from the thermodynamic limit. Obviously in the T-limit different particle moments introduced in Eq. ©, as well as 
the particle number, the total charge and their average values appearing in Eqs. flal) and © are all infinite. Thus, 
to take the thermodynamic limit in Eqs. rial) and © one first expresses the variables (N, Q, N k ) by means of the 
corresponding densities ( n,q,n k ) and then takes the limit V —> oo keeping the densities fixed. This also requires 
the replacement of the discrete sums (I/IO^Cat — > / dn and (1 /V)Y^q f dp by the corresponding integrals over 
densities. 

To formulate correctly the thermodynamic limit of quantities involving densities, one defines the following proba¬ 
bilities 


such that in the limit V —> oo 


P C q {n,V) := VV$ q (Vn,V), 

(10a) 

P f q f(n,q,V) := V 2 V$° q) (Vn, Vq, V), 

(10b) 

Vf q f{q,V) := vv$? q) (Vq,,V). 

(10c) 

P C q {n, V ) = P~(n) + yR c q {n) + 0( V ~ 2 ), 

(lla) 

P f g f(n,q,V) = V%(n,q) + ±Rf q f(n, q )+0(y- 2 ), 

(lib) 

P (9, V) = (q) + i Sf q f(q ) + 0(U- 2 ). 

(11c) 
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The first terms V^ 2 (n) : Wfn) and 'P, c ^!(g) are the T-limits distributions corresponding to the E —> oo limit, the 
second terms in Eq. CD are the finite volume NLO corrections. Obviously, the relation 0 between GC and C 
probabilities also holds in the vicinity of the thermodynamic limit. Thus, the probability distribution ITTbli is just a 
product of C and GC probabilities from Eqs. ill all and lllcl) . 

With the above parametrization of the probability distributions, the particle density moments in the C and the GC 
ensembles are obtained from 

( nk )q = J dnn k/ P^°(n) + J dnn k (n) + 0(V ~ 2 ), (12a) 

( nk )f q f = J dnn k J dqV^(n,q) + ^ J dnn k J dq Rf q f (n,q) + 0{V ~ 2 ). (12b) 

The equivalence of the GC and C ensembles in the thermodynamic limit requires that the first terms in the above 
equations coincide. Indeed it was shown m that the charged density distribution function P/E (g, V ) converges to 
the Dirac delta function such that 


P^)(n,q) =V™(n) ■ S(q- (q)). (13) 

The above relation establishes the equivalence of the GC and C ensembles in the thermodynamic limit on the most 
general, probability level. Indeed, substituting Eq. Cl into Eq. m it is clear that any charged particle density 
moment converges to the same asymptotic value (n k ) oo in the GC and C ensemble. In the T-limit the charge density 
q is also identified with its thermal average value (q) following Eq. 11311 . 

The coincidence of C and GC probability distributions and the corresponding charged particle density moments in 
the asymptotic limit of V —> oo is not any more valid for the NLO contributions. The correction coefficients of order 
(1/E) in Eq. Cl and o are in general different in the C and GC ensembles. 

For physical applications it is important to know the NLO behavior of the probability functions. The asymptotic 
properties of e.g. the semi-intensive quantities will be shown in the next section to be sensitive to the finite volume 
corrections appearing in Eqs. C1C1- 

In the following we will discuss how to obtain the NLO contributions to the probability distributions. We will then 
apply these results to establish the properties of particle density moments as well as semi intensive quantities in their 
approach towards the thermodynamic limit. 


A. Finite volume corrections to the probability distributions and density moments in the C ensemble 

To obtain the finite volume corrections to the probability distribution CD and the corresponding charged particle 
moments CD it is convenient to introduce a generating function 


S(A,E) = ^ A n P$(N,V). 


N=0 


The density moments ( N k ) c are obtained from the generating function (11411 as 


(14) 


(NT = ( A— ) G(X,V) 


A=1 


For the probability distribution in the canonical ensemble 0 the generating function has the following form 

I v q (2Vz 0 VX) 

S(AT) - A iv,(2Vzo) ■ 


(15) 


(16) 


The T-limit behavior of particle moments and the probability distribution is obtained from Eqs. CD and CD using 
the asymptotic behavior of the generating function for V —* oo with fixed q. This is determined by the limiting 
properties of the Bessel function 


lim I a (az) = 

a . —kx) 


y.y/1 -\-z 


V^TTO/I + 2 2 ) 1 / 4 


i + 


i + 


3z 2 - 2 


24(1 + z 2 ) 3 / 2 a 


+ 0(a~ 2 ) 


(17) 
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The generating function m in the T-limit is now obtained as 


Q{\V) = 

V 


q + x q 

Vq | 

_q + x q .\_ 

X 1 


6^A - q 2 6 zq - q 2 „, T __ 2 , 




IZVxl 


+ 0(V~ 2 ) ^ . 


where we have introduced 


= sJq 2 + Az 2 , x q -\ = yjq 2 + Az 2 \ . 


(18) 


(19) 


Applying the above asymptotic form of the generating function in Eq. one establishes (for details see Appendix 
0) the T-limit of positively and negatively charged particle density moments in the canonical ensemble 


k zi 


\k -1 , k{k 1 ) Zq k _ 2 


Af = (n±A - A^{n±) k ~ X + ^W ^^ 2 + 0(V~ 2 ) 


Vx 2 


2V Xn 


( 20 ) 


with x q defined as in Eq. m- The average charged particle densities 


(^±)oo 


Vq 2 +Az 2 ±q 
2 


( 21 ) 


are the asymptotic results obtained from the V —> oo limit and have common values in the C and the GC ensembles 
if one identifies, following Eq. m- the charge density q with its thermal average value (q). 

Taking into account Eq. 11 2 all together with the (1/E) expansion of the charged particle moments 1201 one finds 
that the probability density Pq' (n; V) of negatively charged particles in the near vicinity to the thermodynamic limit 
reads 

P?(n; V) = S (n — {n) ao ) + y^S'(n^ <«>«,) + AA-S"( n - ( ri)oo ) + 0 (V~ 2 ) , (22) 

where the prime and double prime in the delta functions denote the first and the second order derivatives with respect 
to the particle density n. 

Eqs. EH and PH summarize the T-limit behavior of the charged particle density distribution and the density 
moments in the canonical system with exact charge conservation. 


B. Finite volume corrections to the probability distributions and density moments in the GC ensemble 


To establish the asymptotic behavior of the GC probability distribution of charged particle density and the density 
moments in the near vicinity to thermodynamic limit we use the following decomposition of the probability density 
function P^ ( q , n; V) fl^] 


Pffi{q,n-,V) =Pf q f(q;V)pC(n;V). (23) 

The T-limit is obtained from the corresponding behavior of P f^(q\V) and P ^(n\V) distributions. From Eq. 1111) 
one gets 

Pfqf(q,n-,V) = P% ) (q)P?(n) + ± (sgf ( 9 )P|» + R c q {n)V^ («)) + 0(V~ 2 ). (24) 

The coefficient in the (l/E)-expansion of the canonical probability density P^ (n; V) up to 0(V~ 2 ) is already known 
from Eq. PH- However, the GC probability distribution P) 7 !/ (q; V ) is only known in the leading order jl^ | 

Pf g C(q-,V)=S(q-(q))+0(l/V). (25) 

The explicit derivation of the (l/H)-corrections to P^ is given in Appendix B where it is shown that 

Pf,f (?, V) = 8(q-(q)) + ^- 6" (q - (q)) + 0(V ~ 2 ) , 


(26) 
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with X( q ) as in Eq. (ITbll but with q being replaced by (q). 

Substituting the asymptotic expansions ( 1221 ) and ( 12 ( 51 ) to Eq. (1241 one gets after some functional algebra the grand 
canonical probability density as 

p fq) (q, n; V) = 5 {n - {n) 00 )5{q - (q)) + {n - (n) S(q - (q)) + 0 (V~ 2 ) . (27) 

where (n)00 is the average density in the V — » 00 limit which is common for both ensembles. Carrying out the charge 
integration in the above equation one gets the GC particle number density probability distribution 


P'&?(«, V) = 6 (n - <«>«,) + { ^5" ( n - <«)«,) + 0(V ~ 2 ). (28) 

The above result can be directly compared with the corresponding probability distribution Q in the C ensemble. 
It is clear that the leading terms in both ensembles coincide. This is to be expected due to equivalence of the C and 
GC ensembles in this limit. However, the (l/V)-corrections in both ensembles are obviously different. This indicates 
that the C and GC probability distributions converge to the asymptotic, V —> 00 limit with a different strength. 
That is why, in any finite volume the thermodynamic observables calculated in both these ensembles have, in general, 
different values. 

The properties of negatively charged particle number density distributions P^'(n; V) and P;^‘(n, V) calculated by 
means of exact formulae © © are illustrated in Fig© for different volume parameters. It is seen in Fig© that the 
largest differences between C and GC distributions appear for small volumes. With increasing V both distributions 
become narrower, however P'^f' is always broader than . In the large V — > 00 limit both these distributions 
converge to generalized functions which are explicitly described through Eqs. ( 1221 ) and m- 





FIG. 1: Probability distributions of charged particle density in the canonical and grand canonical ensembles calculated for 
different volume parameters (V=l, 10 and 20 in arbitrary units). 

Knowing the asymptotic properties of the probability distributions on can study the corresponding behavior of the 
positively and negatively charged particle density moments. The asymptotic behavior of the positively and negatively 
charged particle density moments in the GC ensemble can be obtained from Eq. ltl2bll and in the following 
transparent form 


<4> GC = (n±) k oo + + 0(V~ 2 ), (29) 

where ( n±} 00 is as before the average density of (±1)- charged particles in the V —> 00 limit which is common for GC 
and C ensemble. In a particular case, for the second moment (fc = 2), the Ea. 12911 multiplied by V 2 is a well known 
result from statistical mechanics: ( N 2 ) — (N) 2 = (N) which tells us that the average square fluctuations in particle 
number in the GC ensemble are controlled by the average number of particles. 


IV. SEMI-INTENSIVE QUANTITIES IN THE THERMODYNAMIC LIMIT 

As an application of the results for the asymptotic behavior of the probability distributions and charged particle 
moments we consider the properties of semi intensive quantities. The simplest example of such a quantity, as discussed 
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in the introduction, is the scaled variance of charged particles 


(Nl) - (iV ±) 2 

- W3 — 


Obviously, one can express u)n ± through the corresponding density moments u n± such as 


(30) 


„(n 2 ± )-(n ± ) 2 
“ N± = -<^>- = Vu ” 


(31) 


The thermodynamic limit of u>n is controlled by the corrections to the first and second moment of the charged 

particle density. From the results in the previous section it is clear that in the T-limit 

(n k ) = (n)^ + 1 U k , q + 0(V ~ 2 ), (32) 

where the first term is common for the GC and C ensembles, however the (1/47) correction term 7 Zk- q differs in these 
ensembles. Indeed, from Eq. (EH and m one finds 


, ''k\q 

V GC 

IK -k-q 


k(k — 1 ) Zq 

2 Xq 


(n) 


k—2 

oo 


k -1 

OO 5 



where TZ C and 1Z GC are the correction terms in C and GC ensemble respectively. 
From Eq. lH21) it is clear that the ratio 


(33) 

(34) 


_ TZk\q 
k ~ V(n)i 


(35) 


describes the NLO corrections to the A:-th order density moment. Following Eqs. m and (El one finds that these 
corrections for the canonical ensemble 


A 


c 

k 


k(k — 1 ) Zq k Zq 

(77)00 I 4 ( 77 )oo 


and for the grand canonical ensemble 


GC _ k(k 1 ) 

^( 77)00 


(36) 


(37) 


are expressed through the total number of particles of a given charge V(n) 00 as well as through dimensionless variables 
Zg/xq^oa and Zq/x 2 which depend only on the q/zo ratio. 

From Eq. GD one finds 


? _ ( rl +)oo , { n -)oo 0 / 00 \ 

2 / \ ' / \ 

2:5 (n_) oo (n+) oo 

thus the NLO corrections A 7 . can be expressed by means of measurable (N_)/(N + ) particles ratios. 

In heavy ion collisions at RHIC charged particle ratios like e.g. 7 r _ / 7 r + , K~/K + or p/p has been measured (see e.g. 

PI ) in the range between 0.6 to 1.0. Thus, for RHIC conditions one estimates from Eq. (1381) that the relevant domain 
of q/zo is between 0 and 0.5. In Fig. E|we calculate from Eqs. (13011 and E3 the NLO corrections to the first, the 
second and the third charge particle density moments for canonical and grand canonical ensemble. These corrections 
are plotted as a function of the multiplicity of negatively charged particles (iV_) at fixed value of (A_)/(7V + ) = 0.8 
particles ratio. 

It is specific for the grand canonical distribution that the first moment Ai coincides with its thermodynamic limit. 
However, as seen in Fig. E this is not anymore the case for higher moments where deviations from thermodynamic 
limit in the GC ensemble are not negligible even for (A_) ~ 150. It is also clear from Fig. Ethat for fixed (iV_) the 
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FIG. 2: Deviations of the first Ai, second A 2 and the third A 3 charged particle moments from the thermodynamic limit as 
functions of charge particle multiplicities (7V_) calculated at fixed (N-)/(N+) = 0.8. The results are shown for the canonical 
and the grand canonical ensembles. 


difference between GC and C ensemble increases with increasing order of particle moments. The convergence of 
to their thermodynamic limit values is slower with increasing k 7 for both ensembles. Finally, it is clear from Fig. [3 
that the first moment is converging from below whereas the second and the third from the above to their asymptotic 
values. 

With the above results for the NLO corrections to different particle density moments one can also study the 
convergence and the thermodynamic limit of the scaled variance PHi. 

Indeed from Eqs. m and 01 one finds: 


T- lim con 


r R.2-,q — 21Zi ;q (n)oo 
(n) 00 


(39) 


Applying in Eq. I® the C and GC correction factors 7 Zk- q from Eqs. PH and one gets in the canonical ensemble 


while in the grand canonical system 


T- lim 1 0 % 


1 yV + 4 zl T q 

2 


T- lim = 1. 


(40a) 


(40b) 


Thus, the scaled variance is finite in the T-limit although it differs in the C and the GC ensembles. This result agrees 
with the previous finding in m and explains how to reconcile this with the T-limit equivalence of different statistical 
ensembles. 

The scaled variance is not the only example of semi-intensive variable. There is actually a broad class of variables 


Sk = 


(. N k ) - ( N) k 


(N) 


k -1 


(41a) 


which have the same properties as ujn'- they are finite in T-limit and have different values dependently on how the 
charge conservation is implemented in the description of the system. Indeed from PH one gets 


T- lim Sk = 


TZk:q klZlqln)^ 1 


(n) 


k -1 


thus, following Eqs. 1331134D the C and GC values for positive(negative) particles in the T-limit are found as 


T-limSf = M*- 1 ) vV+ 4 ^T9 


while in the grand canonical ensemble 


4 vV+4^ ’ 

k(k - 1) 


(41b) 


(42a) 


T- lim Sj? c 


2 


(42b) 
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The scaled variance u>n is just a special case of Sk corresponding to k = 2. 

Another example is a broad class of variables closely related to cumulant or factorial cumulant moments k p defined 
as [ij| 


\p 

In Q(e x , V) = 

V 

\p 

In 5 (A, V) = 


for cumulant moments, 


for factorial cumulant moments. 


In the T-limit those moments are linear in V, thus the ratios 


(43a) 

(43b) 


/C 


p;r — 


Kp 

■> 

K r 


(44) 


are obviously semi-intensive quantities. For cumulant moments the quantity K , 2; i coincides with the scaled variance 
m while for factorial cumulant moments it is equal to lon — 1- 

One can also construct more involved semi-inclusive variables having a finite T-limit behavior which are determined 
by higher order asymptotic terms of the corresponding probability distributions. 


V. SUMMARY AND CONCLUSIONS 

We have considered statistical systems with global charge conservation and discussed the properties of different 
charged particle probability distributions. We have put particular emphasis on the limiting behavior of the probability 
densities in their approach towards the thermodynamic limit. In particular we have calculated the asymptotic values 
and the first subleading finite volume corrections. Our results were obtained in the statistical model of quasi non¬ 
interacting gas of charged particles, constrained by the conservation laws. Such a model has been recently shown to 
be successful in describing particle production in heavy ion collisions and lattice QCD thermodynamics restricted to 
the confined, hadronic phase. 

We have discussed the differences in the asymptotic properties of the probability functions for a system with an exact, 
that is canonical, (C) and with an average, that is grand canonical, (GC) implementation of charge conservation. We 
have shown that in the thermodynamic limit the corresponding probability distributions in the GC and C ensembles 
coincide and are described as generalized functions. This property is a direct consequence of the GC and C ensemble 
equivalence in the thermodynamic limit. However, the first finite volume corrections to the asymptotic value differ 
for both ensembles. 

Finally, using the results of the probability functions we have derived the asymptotic behavior of the charged particle 
moments and established the differences in the GC and C formulation. We have also applied these results to find the 
thermodynamic limit of a class of semi-intensive quantities. It was shown that in systems with exact and average 
charge conservation such quantities should naturally converge to different values in the thermodynamic limit. This is 
because the behavior of the semi-intensive quantities in the near vicinity to the thermodynamic limit are determined 
by the subleading, finite volume, corrections to the probability distributions which are specific to a given statistical 
ensemble. 

Similar conclusions related to the thermodynamic limit in different ensembles can be also drawn for nonrelativistic 
systems H71 . Such an analysis could be of interest for low energy processes in which the particle number conservation 
is preserved. 
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APPENDIX A: LEADING AND NLO COEFFICIENTS OF <1V% 


The leading and next to leading order contributions to the particle moments ( N k ) c are obtained as the coefficients 
Cfc and Cfc_i at large volume expansion terms in Ecp (113, i.e. from 


(N k ) c = Ck + Cfc _ i v fe - i + o(v k ~ 2 ). 


(Al) 


To get these coefficients using the expansion for the Bessel functions given in we first observe that the second 
and the third terms in the curly bracket in Eq. m can be neglected. This is because the contributions of these 
terms to the 0(P fc_1 ) order coefficient exactly cancel each other since the derivative in Eq. 1151) is taken at A = 1. 
Thus, it is sufficient to consider 


Q(X,V) ~ f(X)e v °^ 


(A2) 


with 


/(A) = 


Xq 
%q; X 


g{ A) = x q .x - x q + gin 


q + X q; \ 


where x q is as in Eq. (Hill) . 

To calculate Ck and Cfc_i we use the identity 


dX 


G = Xk d^ g 


k{k 1) k—i d r 
~^ X dX^ G + 


(A3) 


and keep the first two terms with the highest order derivatives. Only the first term in Eq. (IA3I) contributes to Ck 


Ck = f{x) (jx) keVaW • (A4) 

However, the Ck-\ coefficient is receiving contributions from both terms in Eq. (1X1 . The highest derivative term in 
Eq. (HI is calculated from 


d k 


Vg( A) . 


k 

k - 1 


df d 


k— 1 


d\ dX k ~ l 


^Vg(X) . 


where up to 0(V k 3 ) 


e VgW L e Vg(X) +V k-i dk e Vg(X) 

dX k \dX) 

d_ 

dX 

The coefficients dk obey the following recursion relation 


V 


k-1 


k-1 


^ e v B W + v k ~ 2 d k -i e V9 ^ 


d 

dX 


' d k ~ x 

dX^ 1 


e Vg(X) 


(A5) 


(A6) 


dk-dk-^ + ik-l) 


(dg\ k ~ 2 fg 
\dXj dX 2 ’ 


(A7) 


with the solution 

j _ k(k - 1) ( dg\ k ~ 2 d 2 g 
k 2 ^dA/ dX 2 ' 


(AS) 


Applying the above result for dk coefficients in Eq. (IA6I) and using Eq. (IA3I) one gets the hnal expression for the Ck-i 
coefficient 


k (k- 1 ) 

Ck-1 = -x-/(A) 


k -1 


+ k 


df ( dg 


dX \dX 


fc-i 


+ i/w 


A 

dA 2 ’ 


where the last term is obtained from the (fc — l)-order derivative in Eq. 1X1 . 


(A9) 
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APPENDIX B: GRAND CANONICAL CHARGE DENSITY PROBABILITY DISTRIBUTION IN THE 

T-LIMIT 


Let us consider the grand canonical probability distribution to find the charge density q in a system of volume V 
and average charge density (q). This probability, following Eq. (Ill)cl) , reads 


(g) + X{q) 


T Vq 


2zq 


(Bl) 


Pffi(q,V) = Ve- Vx <°U Vq (2Vzo) 

with x/ q ) as in Eq. (11911 but with q being replaced by (q). 

Our goal is to find the first subleading contribution to ( q , V) in the thermodynamic limit. 

It is rather straightforward to see from Eq. m that the limit V —> oo for fixed q does not exist as a regular 
function 0 However, this limit can be found as a generalized function. To show it, we consider the integral 


(G)= J dqG(q)P? q f(q,V). 


(B2) 


where the probability distribution P^f' is smeared out with a test function G. 

In our previous study [12H it was shown that in the T-limit 

(G) = G((q))+0(V~ 1 ). 

Thus, one can write 

lim PfC(q,V)=6(q-(q)) . 

V —►oo 

We apply the asymptotic form of the I a Bessel function ra in m and then calculate the integral 1B2II using the 
saddle-point method. To obtain the NLO term of this integral one should take into account contributions coming 
from the probability EB in the next to leading order and corresponding contributions from (Oil simultaneously. 
From eb and mm one gets 


(B3) 


(B4) 


Pf g f(g,U) ^U 1 / 2 


0 V{ x (q) x q) 




7TX a 


y-1/2 6z 0 - g 2 


12\Z27tq 


1/2 


(g) + X(q) 

q + x q 

(g) + X(q) 


-iVq 


Vq 


V (x {q) -X i ) 


(B5) 


The last term in Eq. contributes to the order 0(V 1 ) in Pf 9 f(g,U) as 

1 6zq — q 2 


V 12v / 27ra;; 


<*(g~ (g» ■ 


A further contribution of order 0(V - 1 ) comes from the coefficient C\ calculated from Eq. PI) as 


ci = 


G(y) f S(y 0 ) - S(y) \ 


-3/2' 


^(y 2 + 4^)V4 V (y-yo) 2 J 


with 


This gives 


S(q) = x q - X( q ) + gin 


(g) + X{q) 

q + Xq 




12 


b (q) 


Consequently, the leading and the subleading contributions to P^ (q, V ) are obtained as 

P(g> V) = S(q- (q)) + ^ S" (q - (q)) + 0(V~ 2 ) , 
where the double prime in the delta function is the second derivative with respect to the charge density q. 


(B6) 


(B7) 


(B8) 


(B9) 


(BIO) 






















12 


APPENDIX C: WATSON-LAPLACE THEOREM 

Let us consider the Laplace integral 

F( A) = J dxe xs{x) f(x). (Cl) 

i 

An asymptotic expansion of F( A) in the limit A —> oo is given by the classical Watson-Laplace theorem: 

Let I = [a, b] be a finite interval such that 

1. max S(x) is reached only in the single point Xq = a, a < xq < b. 

x£l 

2. f(x),S(x) e C(I). 

3. f(x),S(x) € C°° in the vicinity of xq, and S (xo) 7 ^ 0. 

Then, for A £ S e and A —> 00 , the following expansion holds 

OO 

F[ A] ~ e xs{xo) c k \- k ~ 1/2 , (C2) 

fc =0 

where the coefficients 

(C3) 

X=Xq 

S e is here a segment | argt:| ^ ^ — e < ^ in the complex 2 -plane. 


r(fc + 1 / 2 ) ( d\ 2k ( s(x 0 )-s(x) y k ~ 1/2 

(2 k)\ \dx) n> \ (x-x 0 ) 2 J 
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